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Abstract 

Statistical properties of coherent radiation propagating in a quasi - ID ran- 
dom media is studied in the framework of random matrix theory. Distribution 
functions for the total transmission coefficient and the angular transmission 

coefficient are obtained. 
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The discovery of universal conductance fluctuations ( UCF pa has induced a partial 
shift of the main interest in the studies of electronic properties from the averaged values 
of physical quantities to their variance and then to the whole distribution functions (see 
Altshuler et aM and references therein) . Later it was demonstrated that UCF exist also for 
the propagation of classical waves (e. g. light) through disordered systems.^ In contrast to 
electronic measurements which can measure only the conductance of a system, light exper- 
iments have the advantage of being able to measure the angular and the total transmission 
coefficients for an experimental realization. 

In our previous publication! we analyzed the problem of statistics of radiation using 
diagrammatic techniques.lt was rigorously shown that the distribution function can be rep- 
resented through the contribution of connected diagrams only. This representation allowedto 
develop a perturbation theory; in the frameworks of this theory it was found that only for 
moderate values of the angular transmission coefficient the distribution function is a simple 
exponential, as predicted by Rayleigh statistics. For larger values of intensity, the distri- 
bution function differs drastically from a simple exponentialand it's asymptotical behaviour 
is a stretched exponential decay. Also for the total transmission coefficient the Gaussian 
distribution function was obtained. 

An important step was made by Nieuwenhuizen and van Rossum.§ While in Ref. ||]the 
perturbation series was truncated after the second term, Nieuwenhuizen and van Rossum 
using diagrammatic techniques combined with random matrix theory managed to sum up the 
whole perturbation series, obtaining in particular different stretched exponent for the angular 
transmission coefficient distribution function and deviations from the simple Gaussianfor the 
total transmission coefficient distribution function. 

In this paper we reproduce the results of Ref. [5] in the framework of the random matrix 
theory. The approach is based on the analysis of the transfer matrix R (see Stone et aM and 
references therein). Under the restrictions of flux conservation and time-reversal invariance, 
this matrix can be represented in the form: 
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where u and v are arbitrary N x N unitary matrices and A is a real, diagonal matrix with 
N positive elements Ai, . . . , Ajy, where N = W 2 k 2 is the number of transverse channels (W 2 
is the area of the sample). The N x N transmission matrix is given by 

t = u t 1/2 v, (2) 

where r = (1 + A)~ . 

In the isotropic approximation! an ensemble of R matrices is described by the differential 
probability dP(R) = P ({T})Y\dT a dfi(u)dfi(v), where dfi(u) (d/j,(v)) is the invariant mea- 

a 

sure of the unitary group U(N). This isotropic approximation is rather strong assumption 
implying the perfect mode mixing but for a quasi-lD systems it is known to be good.@ 

The angular transmission coefficient T a b, defined as the ratio of the energy carried away 
by the transmitted wave with the transverse wave vector % to the energy of the incident 
wave with the transverse wave vector q a , is given by \t a b\ 2 , t a b being the ah matrix element 
of Eq. (0) . The n-th moment of T a b can be written down in the following way: 

{ T ab) = £ (Km • • • u aan ){u a p x . . . u a p n )*) {(v bai . . . v ban )(v bf3l . . . v bl3n y)o 

M,{/3} 

x{{T ai ...T ari Tp 1 ...Tp n ) 1 ' 2 ) T - (3) 

where the average indicated by the index is performed with the invariant measure of the 
unitary group and (X) T = J d{r} P ({r}) X. 

It is known that to leading order in 1/N both real and imaginary components of u aa and 
Vp b are independently distributed Gaussian variables with zero mean and variance l/2iV.!ll2l 
Then we can write down correlator {{v bai . . . f&an)( u &/3i • • • %3„)*)o as the product ofcorrela- 
tors (v ba v b g) = 5 a p/N summed up with respect to all n\ possible pairings between {a} and 
{/?}. So from Eq. (|3|) we get 
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where 

T a = \Uaa\ T a- (5) 
a 

It can be easily seen, that T a is just the total transmission coefficient: T a = J2bT a b- In fact 
the n-th moment of the total transmission coefficient J2b T a b is 

(\^2 T ab) )= ((u aai ---U aa J(u a / 3l ...U a pJ*) ((v hai ...V bnOln )(v bl p l ...V bn p n )*) 

V b ) W,{/3},w 

x((r ai ...T an T Pl ...Tp n ) 1/2 ) T . (6) 

To leading order in 1/N 

J2((v biai . . . V bnan )(v blPl . . . V bnPn )*) = 8 aiPl . . . 5 anf3n (7) 
{b} 

(because the b-indexes are different we should take into account only one pairing), and the 
right hand part of Eq. (Q) is exactly (T™) . 

Returning to Eq. ([5]) we see that the distribution function P(T a ) can be written as an 
integration over eigenvalues and eigenvectors: 



P(T a ) = fdn.-.j dr n JdU P ({r}) 5 (r a - £ K Q | 2 r^ . 



(8) 



It is convenient to work with the Laplace transform (we also measure T a in units of 
(T a ) = g/N, where g is classical conductance): 

-^exp(sT a ) F{s/g), (9) 

-ioo Z7T2 

Then easilycarrying out the integration with respect to dU, we find: 

N 1 

n*) = ( n t-^t (10) 

a=l a 

We are going to use an approximation of uniform distribution of the "charges" v ai which 
are defined by the relation: r a — 1/ cosh 2 (z/Q,/2).i That is, knowing that the distribution 
of " charges" is statistically homogeneous,! instead of averaging with respect to all possible 
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configurations of "charges" we take into account only one configuration - crystal lattice, 
which leads to the following relation:! 

Ef(r a )=gf \ % f(r). (11) 
±^ Jo 2rVl - r 



for any /(r) which goes to zero when r goes to zero. Then fromEq. fliTf ) we get 



F(s) = exp 



-g 



f \ 1 ln(l + ar) = exp\-g\n 2 (VTT~s + ^)] , (12) 
Jo 2rv 1 — t - 



which exactly coincides with the result of Ref. || Eq. (0) gives in particular Gaussian 
behavior for T a sa 1: 

P(Tj^^exp[-^(T a -l) 2 ] (13) 

and simple exponential decay for large T a : 

P(T a ) ~ exp(-#T a ). (14) 

Now let us return to Eq. (^). As it is known! it means 

P{T ab ) = | o °° dT a P{T a ) i exp (- (15) 

(we measure T afe in units of (T ab ) = g/N 2 , where g is classical conductance). This distribution 
functioncan be described as the Rayleigh distribution function for the angular transmission 
coefficient but with some effective averaged value which inturn fluctuates around the real 
averaged value and the latter fluctuations are described by the total transmission coefficient 
distribution function.!'! Eq. (|i"5|) gives in particular Rayleigh statistics 

P(T ab ) « exp (-T ab ) . (16) 

for T ab <C sjg, ! and stretched exponential tail 



P(T ab )~exp[-2^gT ab ). (17) 

for > 9$ 



Having in mind the comparison of the theoretical result with an experiment it is con- 
venient to express g through the first two moments either of the totalor of the angular 
transmission coefficient distribution function. Calculating the coefficient before s 2 in the 
expansion of the exponent in the right hand part of the Eq. (1T2|) we get: 



< T 2 > 2 

1 = TT- , (18) 



< T a > 2 3g 
and 

< T ab > 2 3g 

which exactly coincides with the result of Ref. §. 

In conclusion we want to discuss the difference between the statistics of total transmission 
coefficient and statistics of conductance g = J2m T m- Taking into account the bimodal 
distribution of r we may say, at least qualitatively, that the conductance is simply the number 
of " open" channels05 g = N e ff . The total transmission coefficient is also the sum with 
respect to "open" channels but each channel comes with a random weight. So the Gaussian 
law for the total transmission coefficient distribution function is just the manifestation of 
the Central Limit Theorem, which is true when N e ff — > oo. In the paper we are taking into 
account the finitness of theparameter N e ff, which is important in particular for obtaining 
correct asimptotics. On the other hand, the conductance fluctuations are determined by the 
strongly suppressed fluctuations of the number of open channels, which in our case can be 
neglected. This principal difference between the two statistics wouldalso manifest itself if one 
tries to go beyond quasi- ID. While the eigenvalue distribution (and hence the conductance 
distribution function) can be not very sensitive to the dimensionality and stay bimodal as 
long as we are in a diffusive regime,0 the isotropic approximation which was essential in 
obtainingEq. flTj) ceases to be valid beyond quasi- 1D.0 
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